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INTRODUCTION 


Problem  Statement 

This  fellowship  addresses  two  important  steps  in  the  battle  against  breast  cancer:  earlier 
detection  through  improvements  in  diagnostic  techniques  and  improved  quality  control  through 
training  of  qualified  personnel,  specifically  diagnostic  radiology  physicists.  The  severity  of  the 
impact  of  mammary  carcinoma  on  the  health  and  mortality  of  US  women  has  been  described 
extensively  in  the  literature.  Many  approaches  to  the  improvement  in  the  overall  treatment  of 
the  problem  have  been  proposed  and  are  being  pursued.  Among  approaches  to  address  this 
problem  is  research  in  improved  breast  imaging  techniques  for  earlier  detection.  The  latter 
stresses  the  need  for  the  general  availability  of  safe  and  low-cost  diagnostic  techniques  to 
monitor  the  population  at  risk.  This  is  particularly  urgent  in  view  of  the  fact  that  the  persistent 

mortality  rate  despite  improvements  in  diagnosis  and  treatment  is  in  part  explained  by  the  fact 
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that  minority,  older  and  low-income  women  have  restricted  access  to  mammography.  While 
conventional  mammography  is  quite  effective  as  a  screening  tool  for  the  detection  of 
microcalcifications,  early  detection  of  breast  cancer  is  frequently  complicated  by  the  occurrence 
of  benign  fibrocystic  changes.^’®  These  changes  include  cyst  formation,  calcium  and  cholesterol 
precipitation,  lobular  secretion  and  epithelial  hyperplasia  and  are  normal  variations  of  breast 
tissue  growth.  The  formation  of  these  benign  lesions  most  commonly  occurs  in  the  fatty  tissue 
which  surrounds  the  glandular  tissue  and  these  lesions  generally  have  smooth,  circumscribed 
borders.  This  normal  dynamic  evolution  of  the  breast  increases  the  probability  of  incorrect 
diagnosis  of  early  mammary  carcinoma. 


Bacl^round 

The  accreditation  of  medical  physicists  in  Diagnostic  Imaging  Physics  and 
Mammography  by  professional  organizations  such  as  the  American  Association  of  Physicists  in 
Medicine  (AAPM)^,  the  American  College  of  Radiology  (ACR)^,  the  Federal  Drug 
Administration  (FDA)  or  the  State  of  California  (Dept,  of  Health  Services)  has  general  education 
as  well  as  experience  requirements.  The  former  were  addressed  through  graduate  studies  in  the 
Department  of  Radiological  Sciences  at  the  University  of  California,  Irvine,  and  through 
participation  in  specific  professional  training  courses.  Work  experience  requirements  were 
principally  met  through  my  work  as  medical  physicist  in  the  University  of  California,  San  Diego, 
Department  of  Radiology  and  my  woe  (work  without  compensation)  association  with  the 
Veterans  Affairs  Medical  Center  in  La  Jolla,  CA. 

The  standard  for  diagnostic  breast  imaging  has  evolved  over  many  years  to  be  X-ray 
mammography.  It  has  superb  resolution  and  clarity  for  detection  of  clusters  of  microcalcifications 
and  for  identification  of  lesion  borders.  However,  early  carcinoma  is  often  difficult  to  distinguish 

from  benign  fibrocystic  changes  and  this  leads  to  a  high  number  of  false  positive  diagnoses  and 
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thus  unnecessary  biopsies.  ’  There  has  been  extensive  research  to  provide  adjunctive  or 
alternative  imaging  systems  to  mammography  including  digital  x-ray  imaging,  magnetic  resonance 
imaging  (MRI)  and  ultrasound  sonography.  In  the  early  1980’s,  dedicated  ultrasound  scanners  were 
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evaluated  for  whole  breast  imaging.  .  Today,  the  role  of  ultrasound  in  breast  imaging  is 
primarily  limited  to  evaluating  the  composition  of  larger  lesions  detected  by  mammography  in 

terms  of  cystic  (benign)  versus  solid  (benign  or  malignant).^  ^ 

Computed  tomography  (CT)  is  a  powerful  and  routine  technique  in  x-ray,  nuclear  medicine 
and  magnetic  resonance  imaging.  In  the  1970’s  and  early  1980’s,  several  investigators  sought  with 

limited  success  to  apply  these  image  reconstruction  methods  to  pulse-transmission  ultrasound  of 
12-17 

breast  tissue.  The  goal  of  these  methods  was  to  improve  differentiation  of  breast  cancer  from 
numerous  benign  conditions  of  the  breast  by  measuring  several  acoustic  properties  of  the  lesion.  A 
particularly  interesting  approach  to  ultrasoimd  CT  is  diffraction  tomography  which  has  received 
substantial  theoretical  treatment  and  holds  the  potential  to  account  for  the  inherent  diffraction  in 

sound  propagation.  Despite  considerable  theoretical  development  of  diffraction  tomography, 
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laboratory  implementation  has  been  limited  and  has  met  with  variable  success. 

The  work  which  constitutes  the  dissertation  research  portion  of  this  fellowship  is  based  on 

a  new  approach  to  ultrasound  diffraction  tomography  which  is  largely  the  result  of  a  collaboration 
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between  researchers  at  UCSD  and  Thermotrex  Technologies  Corporation  (TTC).  This 
approach  is  described  briefly  in  a  subsequent  section  (Methods). 


Purpose 

Certification  of  medical  physicists  in  Diagnostic  Radiological  Physics  and  Mammography 
permits  them  to  assure  and  verify  the  proper  and  safe  performance  of  diagnostic  imaging 
equipment  and  thus  protects  the  safety  of  the  patients.  Additionally,  maintenance  of  imaging 
performance  standards  results  in  improved  and  more  accurate  diagnostics  and  thus  serves  the 
general  patient  population.  This  is  particularly  critical  in  mammography  where  the  imaging 
procedures  have  to  be  optimized  to  permit  earliest  detection  of  breast  cancers.  This  has  recently 
been  recognized  by  the  establishment  of  more  stringent  standards  both  for  the  certification  of 
medical  physicists  and  the  performance  of  mammography  centers.  Thorough  training  in  diagnostic 
imaging  physics  and  mammography  for  medical  physicists  thus  proves  to  be  a  benefit  to  the  local 
community  in  that  the  skills  and  training  acquired  will  serve  the  general  population  at  risk,  as  well 
as  older,  minority  and  poor  women  who  presently  receive  inadequate  diagnostic  screening  for 
breast  cancer. 

The  principal  purpose  in  the  development  of  the  ultrasound  CT  scanner  by  the  UCSD/TTC 
team  was  that  the  differences  described  above  between  this  new  approach  and  previous  efforts  to 
construct  a  practical  ultrasound  CT  (USCT)  system  might  prove  to  be  advantageous.  The  prineipal 
advantages  were:  1)  Low-frequency,  continuous- wave  insonification;  2)  large  number  of 
projections  around  the  whole  object;  3)  high  speed  sampling  at  optimal  spatial  frequency;  4)  large 
field  of  view.  A  laboratory  prototype,  described  in  the  Methods  section,  suggested  that  the 
approach  may  have  utility  as  an  adjunctive  tool  to  mammography.  The  human-sized  device  which 
was  developed  produces  unique  images  of  a  very  large  eight-inch  field  of  view  in  a  coronal 
format.  In  cases  where  an  x-ray  mammogram  reveals  a  suspicious  lesion,  monitoring  with  USCT 
could  provide  a  non-invasive  alternative  to  biopsy. 

The  focus  of  the  dissertation  research  is  the  quantitative  description  of  the  imaging 
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characteristics  of  this  new  ultrasound  computed  tomography  procedure.  Diffraction  tomography 
potentially  affords  the  characterization  of  certain  breast  tissue  properties  such  as  speed  of  sound, 
density,  attenuation  and  compressibility.  Before  a  quantitative  interpretation  of  such  parameters 
can  be  realized,  however,  the  system  needs  to  be  calibrated  with  phantoms  that  establish  the 
linearity,  resolution  and  dynamic  range  of  the  reconstructed  images.  This  is  particularly  essential 
in  the  development  of  a  clinical  breast  imaging  system  since  breast  tissue  is  complex,  the  breast 
has  few  recognizable  anatomical  landmarks  and  thus  imaging  artifacts,  noise  and  aberrations 
must  be  carefully  characterized  and,  if  possible,  eliminated,  in  order  to  reduce  the  probability  of 
mis-diagnosis  or  misinterpretation  of  the  images. 


Methods 

The  emphasis  of  this  pre-doctoral  fellowship  is  on  two  areas  that  are  related  in  subject 
matter  but  rather  distinct  in  terms  of  the  activities  and  approaches  involved.  The  first  task 
addresses  preparation  for  board  certification  in  Diagnostic  Radiological  Physics  and 
Mammography  and  thus  contains  both  formal  and  on-the-job  training  elements.  Certification  in 
diagnostic  radiological  physics  requires  both  formal  training  leading  to  a  graduate  degree  in 
physics,  medical  physics  or  a  related  subject  area  as  well  as  documented  work-experience  in  the 
field  under  the  supervision  of  certified  medical  physicists.  This  report  describes  the  training, 
work  experience  and  degree(s)  conferred  during  the  first  year  of  this  fellowship.  The  second  task 
involves  dissertation  research  on  a  fiill-breast  imaging  technique  based  on  ultrasound  diffraction 
tomography  principles  realized  in  a  recently  developed  scanning  system. 

1.  Certification/Diagnostic  Radiological  Physics  and  Mammography 

The  AAPM  Report  on  the  Clinical  Training  of  Radiological  Physicists’  defines  the  areas 
of  activity  of  the  Diagnostic  Radiological  Physicist  as  follows:  1)  Calibration  of  imaging 
equipment;  2)  Calculation  and  measurement  of  exposure  and  dose;  3)  Improving  and 
maintaining  medical  image  quality;  4)  Training  of  physicists,  clinical  diagnostic  imaging 
residents,  radiological  and  ultrasound  technologists,  and  other  allied  health  professionals  in 
diagnostic  radiology;  5)  Education  of  health  professionals  and  the  public  in  diagnostic  imaging 
physics  and  radiation  effects;  6)  Additional  duties  such  as  developmental  studies,  administrative 
duties  and  research. 

Competence  in  the  above  radiological  and  clinical  areas  derives  from  fundamental 
understanding  of  and  facility  with  radiation  physics,  general  and  specific  imaging  concepts, 
radiation  biology  and  basic  clinical  concepts  including  anatomy,  physiology  and  patient 
procedures. 

Effective  January  1,  1993,  The  American  College  of  Radiology  (ACR)  has  established 
qualification  requirements  for  accreditation  of  medical  physicists  in  mammography  under  the 
Mammography  Quality  Standards  Act  (MQSA).  These  specify  that  for  eligibility,  one  of  two 
sets  of  criteria  must  be  met: 

1)  Medical  physicist  certified  in  Diagnostic  Radiological  Physics  or  Radiological  Physics 
by  any  certifying  body  recognized  by  the  American  College  of  Radiology. 
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2)  Until  January  I,  1996,  alternative  criteria  will  be  recognized. 

a)  M.S.  or  M.A.  degree  in  physics,  applied  physics,  radiological  physics, 
biophysics,  health  physics  engineering,  public  health  when  B.S.  is  in  a  physical  science. 

b)  Training  in  biological  sciences. 

c)  Minimum  1  year  training  in  medical  physics  in  the  area  of  diagnostic 
radiological  physics. 

d)  Minimum  2  years  experience  conducting  mammography  equipment 
performance  evaluations. 

Additionally,  it  is  specified  that  a  candidate  satisfying  either  of  the  above  requirements 
have  had  at  least  fifteen  (15)  hours  of  documented  continuing  medical  education  training  in 
mammography  physics  within  the  last  three  years. 

2.  Dissertation  Research 

The  focus  of  the  doctoral  dissertation  is  the  quantitative  description  of  the  imaging 
characteristics  of  a  new  ultrasound  computed  tomography  procedure  that  is  intended  to  serve  as 
an  adjunct  to  conventional  x-ray  mammography.  This  new  approach  is  based  on  two  human-size 
systems  designed  specifically  for  breast  imaging  that  use  low-frequency  (0.2-1 .4  MHz) 
continuous- wave  sound ,  cylindrical  arrays  of  1024  and  512  transducers,  high  spatial  frequency 
sampling  of  the  wavefront  and  a  unique  diffraction  tomographic  reconstruction  method. 

This  method  of  ultrasound  CT  employs  a  full  circle  of  transducers  which  surrounds  and 
defines  the  field  of  view  (Figure  1).  The  object  is  coupled  to  the  transducer  array  by  a  water 
bath.  One  element  at  a  time  transmits  a  continuous  single-frequency  sound  beam  with  a 
cylindrical  wavefront.  The  system  reaches  steady  state  in  about  300  psec,  at  which  point  the 
remaining  transducers  record  amplitude  and  phase  (relative  to  the  transmitter)  of  the  signal 
scattered  the  object  to  be  imaged .  Each  transducer  is  activated  as  a  transmitter  in  turn, 
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sequentially  probing  along  a  complete  360  arc. 

The  laboratory  systems  were  constructed  to  test  feasibility  of  the  method  for  breast 
imaging.  The  1  MHz  SCT  device  shown  in  Figure  2  consists  of  1024  lead  zirconate  titanate 
(PZT)  transducers,  arranged  in  an  eight-inch  cylinder.  The  breast  imaging  system,  shown  in 
Figure  3,  consists  of  a  control  console  which  contains  all  the  signal  generation  and  data 
acquisition  and  processing  electronics  and  computers  and  a  couch  on  which  the  patient  reclines 
with  the  breast  suspended  in  a  water  tank  incorporating  the  transducer  array.  The  tank  can  be 
raised  and  lowered  for  slice  selection  and  moved  laterally  for  optimum  patient  positioning. 

The  emphasis  of  the  first  year’s  work  was  on  evaluating  the  acoustical  characteristics  of 
the  two  transducer  rings,  on  establishing  initial  system  performance  figures  of  merit,  on 
assessing  the  effects  of  various  object  and  data  acquisition  variables  on  image  quality  and  on 
participating  in  a  clinical  trials  series.  These  tasks  followed  the  plan  in  the  Statement  of  Work  of 
the  original  proposal.  An  additional  task  was  introduced  in  the  first  year  and  expanded  in  the 
second:  analysis  of  the  raw  signals  scattered  by  breast  tissue.  This  work  is  based  on  similar 
statistical  analyses  reported  by  Steinberg  and  co-workers  as  well  as  other  investigators^^'^^  who 
have  studied  the  amplitude  and  phase  distortions  of  the  wavefront  of  the  insonifying  beam  when 
passing  through  tissue  and  frequently  breast  tissue.  The  relevance  of  this  sort  of  data  analysis  is 
that  it  describes  and  summarizes  the  strength  of  the  scattering  of  the  ultrasound  beam  and  thus 
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was  felt  to  be  potentially  an  indicator  of  the  suitability  or  lack  thereof  of  certain  reconstruction 
algorithms  based  on  linearizing  assumptions  such  as  the  Bom  approximation.  A  brief  discussion 
is  presented  of  a  new  operator-based  formalism  that  suggests  that  the  Bom  and  Rytov 
approximations  represent  extreme  endpoints  of  a  continuos  reconstmction  formalism.  Due  to 
limited  access  to  the  imaging  system  as  a  result  of  equipment  problems  in  the  second  year,  a 
second  new  task  was  added,  the  analysis  of  breast  tissue  contrast  for  patient  data  available  from 
clinical  trials  conducted  at  the  UCSD  Center  for  Women’s  Health.  Both  normal  tissue  and  tissue 
diagnosed  with  pathology  on  the  basis  of  conventional  mammography  were  evaluated. 

The  next  section  contains  a  description  of  the  work  performed  during  the  fellowship  and 
examples  of  some  of  the  results  obtained,  as  well  as  discussions  of  the  experimental  and  analysis 
techniques  employed.  The  tasks  addressed  are  summarized  below: 

•  Characterization  of  the  transmit  beam  patterns  and  slice  profiles  of  both  the  1 
MHz  and  0.5  MHz  SCT  rings. 

•  MTF  (Modulation  Transfer  Function)  measurements  vs.  frequency  and  position 
within  the  measurement  tank  for  the  1  MHz  SCT  ring. 

•  Speed  of  Sound  contrast  linearity  for  various  size  objects  in  a  non-aberrating 
medium  (water)  for  both  1  MHz  and  0.5  MHz  SCT  rings  determined  with  several  reconstruction 
algorithms. 


•  Studies  of  in  vivo  scattering  amplitude  distribution  of  breast  tissue  from  patients 
scarmed  during  the  clinical  trials  phase. 

•  Exploration  of  a  new  operator-based  formalism  for  the  Bom  and  Rytov 
approximations. 

•  Breast  tissue  contrast  and  texture  analysis  from  clinical  patient  data  obtained  at 
the  UCSD  Center  for  Women’s  Health. 

•  Presentation  of  papers  at  the  AAPM  Annual  Meetings  in  Anaheim,  CA  in  July, 
1994,  and  Philadelphia  in  July,  1996;  at  the  American  Institute  of  Ultrasound  in  Medicine 
Conference  in  San  Francisco,  CA  in  March,  1995,  and  at  the  22nd  International  Acoustical 
Imaging  Conference,  Florence,  Italy,  Sept.  1995. 
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Figure  1.  Scanning  Geometry  of  Ultrasound  Diffraction  Tomography  Breast 
Imaging  System. 
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♦  1024  transducers,  20  cm  bore 

•  Operating  frequency, 700  kHz  -  1 .2  MHz 


•  0.37  mm  image  resolution 


•  Data  acquisition  time,  0.5  sec 

•  Image  reconstruction  time,  2  minutes 
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Figure  2.  1024-element  1  MHz  SCT  Transducer  Ring. 
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CLINICAL  TEST  UNIT  LO^iAD 


Figure  3.  TTC/LORAD  SCT  Clinical  Test  Unit 
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BODY 


e  mam  accomplishments  of  the  pre-doctoral  fellowship  in  both  the  training  and  dissertation 
research  areas  are  described  in  this  section.  Results  of  the  research  on  the  performance 
evaluations  of  an  ultrasound  diffraction  tomography  breast  imaging  system  and  of  some  early 
c  inical  patient  data  are  summarized  The  thrust  of  this  report  is  not  only  to  present  final  results 
and  conclusions  about  the  technology  issues  raised  in  the  proposal  for  this  work,  but  to  indicate 

progress  made  and  to  point  out  areas  which  are  new  or  which  differ  significantly  from  the 
statement  of  work. 


I. 


Certification  in  Diagnostic  Radiological  Physics,  Mammography 


Below  are  summanzed  the  training  activities,  formal  education  elements  and  practical 
work  expenence  that  all  satisfy  various  requirements  for  board  certification  in  Diagnostic 
a  lological  Physics  and  m  Mammography  by  agencies  such  as  the  American  Association  of 
Physicists  in  Medicine(AAPM),  the  American  College  of  Radiology  (ACR),  the  American 
College  of  Medical  Physics  (ACMP),  the  State  of  California  (Dept,  of  Health  Services)  and  the 
Federal  Drug  Administration  (FDA). 

1.1  Formal  Education 


Course  Work-Three  academic  quarters  of  Clinical  Rotation  in  several  areas  addressing 
vanous  magmg  modalities  were  completed  in  the  Radiological  Sciences  Department  of  the 
mversify  of  California,  Irvine,  CA  (UCI).  These  courses  were  conducted  mostly  at  the  UCI 
Medical  Center  m  Orange,  CA  but  also  included  visits  to  imaging  centers  such  as  the  UCI  Brain 
maging  Center.  These  rotetions  emphasized  observation  of  and  participation  in  daily,  routine 
activities  in  a  hospital  setting  that  related  to  the  operation,  quality  control  and  interpretation  of 
various  imaging  systems.  The  rotations  were  as  follows: 

•  Radiography^luoroscopy-Participation  in  system  quality  control 

activities,  acquisition  of  images,  review  of  x-rays  with  radiologists. 

f  .  *  X-Ray  CT-Cahbration  of  a  CT  Scanner  using  a  Calphan  phantom,  measurements 
of  system  resolution,  contrast  sensitivity,  slice  thickness,  linearity.  Observation  of  patient  scans 
and  review  of  tomographic  images  with  radiologists. 

MRI-Clinical  applications  presentation  by  radiologists,  observation  of  patient 
scans,  review  of  patient  images  with  radiologists,  daily  and  weekly  system  calibration  procedures 
using  several  phantoms. 

Mammography/Ultrasound-Participation  with  sonographers  on  patient  scans 
evaluation  of  mammograms  and  adjunctive  ultrasound  scans  of  a  range  of  breast  types  and 
patholopes  with  radiologists,  daily  ultrasound  and  mammography  systems  checks  and 
calibrations. 
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Radiation  Oncology-Review  if  linac  calibration  procedures,  generation  of 
patient  depth-dose  profiles  from  bremsstrahlung  and  electron  beam  irradiation,  attendance  of 
weekly  tumor  review  boards,  review  of  treatment  planning  procedures  and  participation  in 
patient  treatments. 


•  Nuclear  Medicine-Presentation  of  and  participation  in  radiopharmaceuticals 
preparation  and  assaying,  observation  of  gamma  camera  patient  imaging,  daily  camera 
calibration  procedures,  review  with  radiologists  of  SPECT  images,  review  of  PET  procedures 
including  cyclotron  operation  for  radiopharmaceutical  production,  observation  of  patient 
treatment  and  imaging,  review  of  images  with  radiologists. 

Digital  Angiography— Review  of  digital  fluoroscopy  procedures,  participation  in 
system  performance  calibrations,  implementation  of  various  digital  subtraction  techniques. 


Degrees  Conferred:  M.S.  degree  in  Radiological  Sciences,  March,  1995.  University 
of  California,  Irvine,  CA. 

Doctoral  Dissertation:  “The  Empirical  Determination  of  the  Range  of  Validity  of  the 
Bom  Approximation  in  Ultrasound  Diffraction  Tomography”  with  Prof  Sidney  Leeman 
University  of  London,  Kings  College,  London,  UK. 


1.2  Training/Clinical  Experience 


•  Mammography  Practicum-Participated  in  a  1  week,  40  hour  practicum  for 
mammography  technologists  offered  by  the  UCSD  School  of  Medicine  at  the  UCSD  Medical 
Center.  The  course  consisted  of  both  lectures  and  working  sessions  and  emphasized  breast 

pathology,  optimal  mammography  techniques  and  review  and  interpretation  of  mammograms 
with  radiologists. 

rncT-v  *  Clinical  Trials-Participated  in  the  performance  of  clinical  trials  of  patients  at  the 
UCSD  Center  for  Women's  Health  designed  to  correlate  conventional  mammograms  and 
sonography  with  ultrasound  CT  images.  About  twenty-five  patients  with  a  wide  range  of  breast 
types  and  pathologies  were  studied.  The  work  was  carried  out  in  close  collaboration  with 
mammography  radiologists. 

•  Diagnostic  Imagin^ammography  Quality  Control-Participated  in  a  wide- 
ranging  number  of  diagnostic  imaging  equipment  quality  control  surveys  and  radiation  shielding 
planning  and  evaluation  procedures.  These  activities  are  summarized  in  a  table  in  Appendix  1. 
Of  particular  importance  was  participation  in  the  acceptance  testing  of  a  GE  Senographe  DMR 
Mammography  unit  at  the  Veterans  Affairs  Medical  Center,  La  Jolla,  CA.  in  Dec.  1993  and  the 


9 


subsequent  submission  of  calibration  and  performance  data  to  the  ACR  for  facility 
accreditation.  Accreditation  was  conferred  on  the  first  attempt  in  February,  1995. 

Mammography  Courses/Seminars-Attended  courses  dealing  with 
mammography  quality  control  practices  and  regulations  both  at  the  AAPM  Annual  Meetings  in 
^aheim,  CA  m  July,  1994  and  in  Philadelphia,  PA  in  July,  1996  as  well  as  at  Southern 
Cahiomia  Chapter  AAPM  meeting  in  Carlsbad,  CA,  May,  1995. 


•  AAPM  Summer  School-Received  29  CEU  units  in  Brachytherapy  Physics  by 
attending  the  1 994  AAPM  Summer  School  at  UCSD,  San  Diego,  CA  in  June,  1 994. 


•  Mammography  Quality  Standards  Act  (MQSA)-Intensive 
offered  by  AAPM,  12  MPCEC  units.  Nov.  3-4, 1995,  Pasadena,  CA. 


Practical  Review 


f  I  and  Interventions  Course- 1 5  hours  of  CEU  by  the  University 

ot  California  at  San  Diego  School  of  Medicine  April  12-14,  1996  in  breast  imaging  and 
mte^entions  including  diagnosis  of  breast  pathology,  quality  standards  and  reporting  fine- 
needle  core  biopsies  and  ductography. 

.1,  c.  /  Mammography  Physicist-Submitted  application  package  to 

the  State  of  California  for  registration  as  mammography  physicist.  Part  of  the  application  is  a 
complete  quality  control  survey  of  a  mammography  unit  which  was  a  GE  DMR  unit  at  the  VA 
Medical  Center  in  La  Jolla,  CA. 


II.  DISSERTATION  RESEARCH 

Tasks  consistent  with  the  Statement  of  Work  of  the  fellowship  proposal  are  briefly 
summarized.  Example  results  are  presented  and  some  preliminary  conclusions  are  drawn  where 
appropnate.  One  task  which  was  not  called  out  directly  in  the  proposal  is  the  statistical 
evaluation  of  the  raw  data  scattered  by  breast  tissue.  This  follows  recent  work  reported  in  the 
literature  where  such  analyses  are  used  to  correct  the  transducer  wavefront  distortions  in 
conventional  B-scan  ultrasound  imaging.  The  application  here,  which  extends  the  data  to  much 
lower  frequencies  than  those  reported  before,  is  to  use  this  analysis  to  characterize  breast  tissue 
an  to  0  tain  an  a  priori  indicator  of  the  degree  of  scattering  and  thus  the  validity  of  particular 
reconstruction  algorithms.  A  second  task  which  was  added  to  the  scope  of  the  proposed 
Statement  of  Work  is  the  evaluation  of  breast  tissue  contrast  for  some  of  the  patients  imaged  as 
part  of  the  clinical  trials.  This  task  was  added  when  it  became  clear  that  some  of  the  originally 
proposed  system  performance  characterizations  could  not  be  completed  due  to  problems  with 
systems  components  that  significantly  reduced  the  availability  of  the  SCT  (Sonic  Computed 
Tomography)  breast  imaging  system. 

A  further  expansion  of  the  SOW  is  represented  by  the  fact  that  during  the  reporting  period 
a  second  transducer  ring  centered  at  0.5  MHz  was  fabricated  and  integrated  with  the  system.  The 
motivation  for  adding  such  a  low  frequency  ring  arose  out  of  the  clinical  trials  which  indicated 
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substantial  attenuation  and  phase  shifts  for  either  large  or  very  dense  breasts.  A  lower  fi-equency 
ring,  though  resulting  in  loss  of  resolution,  was  expected  to  have  better  signal  to  noise 
characteristics  and  thus  provide  improved  images.  Many  of  the  tasks  listed  in  the  Statement  of 
Work  have  been  performed  for  both  transducer  rings. 


2.1  Transducer  Beam  Pattern  Characterization 


One  of  the  most  fundamental  system  characterizations  is  the  determination  of  the 
transducer  radiation  pattern  and  the  effective  slice  width.  The  radiation  pattern  must  be  specified 
accmately  for  solution  to  the  wave  equation.  For  the  sake  of  speed  of  calculation  and 
mat  ematical  simplicity,  it  is  presently  assumed  to  be  represented  by  the  first  derivative  of  the 
Hankel  function,  that  is,  it  is  assumed  to  be  dipole  radiation  with  a  cos0  distribution.  The  extent 
o  which  the  actual  radiation  (and  receive  )  pattern  deviates  from  such  a  distribution  is  of 
fundamental  interest. 

Measurements  were  carried  out  with  the  experimental  geometry  sketched  in  Figure  4  A 
needle  hydrophone  with  approximately  0.03  in.  diameter  transducer  element  was  positioned ' 
opposite  an  active  transducer  at  roughly  the  center  of  the  tank.  The  needle  hydrophone  was 
translated  across  the  tank  diameter  and  at  each  measurement  position  it  was  rotated  to  maximize 
signal,  thus  assuring  that  the  position  was  radial.  A  second  measurement  series  was  carried  out 
by  recording  the  on-axis  beam  intensity  as  a  function  of  distance  from  the  transducer.  This 
information  IS  not  only  useful  for  mapping  out  the  near  to  far-field  transition  but  was  used  to 
correct  the  beam  profile  data  to  constant  radial  distance  from  the  transducer.  For  the  0.5  MHz 
transducer  nng,  the  needle  hydrophone  was  found  to  be  too  insensitive  to  record  the  SCT  (sonic 

computed  tomography)  transducer  signals  and  a  larger,  0.24  in.  diameter  transducer  was 
employed. 

Figure  5  is  a  polar  plot  of  the  radiation  pattern  for  one  of  the  transducers  of  the  1  MHz 
SCT  nng.  It  is  compared  with  both  a  cosO  and  a  (cosO  f  distribution  and  is  seen  to  fit 
somewhere  in-between.  The  situation  for  the  0.5  MHz  ring  is  depicted  in  Figure  6,  which  shows 
at,  even  with  the  coarseness  of  the  data,  the  beam  pattern  is  characterized  by  side-lobes  that 
raise  serious  questions  about  the  effect  of  the  dipole-radiation  assumption  on  the  image 
reconstmction  The  results  for  the  lower  frequency  ring  were  anticipated  since  the  lower  length 
to  width  (L^)  ratio  (due  to  greater  detector  width,  1.2  mm  vs.  .6  mm  for  the  1  MHz  ring)  results 
in  a  non-cylmdncal  antenna  pattern. 
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Rel.  AmpHtude 


Figure  4.  Schematic  Diagram  of  SCT  Transducer  Beam  Pattern  and  Slice  Profile 
Determination  using  Needle  Hydrophone 


1  MHZ  RING 


Figure  5.  Polar  Plot  of  Transducer  Radiation  Pattern.  1  MHz  SCT  Ring 
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RADIATION  PATTERN 
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Figure  6.  Radiation  Pattern  for  0.5MHz  SCT  Ring  Showing  Side  Lobes. 


2.2  Tomographic  Slice  Geometry 


The  measurement  of  effective  slice  geometry  was  performed  at  the  center  of  the  tank 
using  the  needle  hydrophone  for  the  1  MHz  ring  and  the  0.25  inch  dia.  transducer  for  the  lower 
frequency  ring.  One  transducer  was  excited  and  the  axial  pulse  intensity  distribution  was 
measured  by  vertically  translating  the  recording  transducer.  The  results  for  the  high  frequency 
ring^e  displayed  in  Figure  7  and  suggest  that  the  slice  width  is  essentially  that  of  the  transducer 
eight  ^),  though  the  profile  for  the  0.5  MHz  ring  is  not  nearly  as  well  defined  as  that  for  the 

higher  frequency  ring.  It  is  to  be  noted  that  these  data  obtain  only  for  one  transducer  at  one 
location,  the  center  of  the  tank. 

What  IS  also  of  interest  is  the  effective  slice  width  or  sensitivity  profile,  measured  as  the 
axial  range  over  which  information  shows  up  in  the  reconstructed  image.  Such  a  measurement 
was  performed  with  two  orthogonal  nylon  filaments  inclined  at  approximately  45”  angles.  The 
filaments  were  120  pm  in  diameter  and  were  too  small  to  be  imaged  with  the  0.5  MHz  SCT  ring 
With  the  1  MHz  ring,  however,  the  resultant  images  were  1 . 1  and  1 .2  cm,  respectively,  very 
much  in  agreement  \vith  the  width  of  the  profile  in  Figure  7. 
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SONIC  CT  SLICE  PROFILE 


Figure  7.  Slice  Profile  for  1  MHz  SCT  Ring 


2.3  System  MTF  Measurements 

K  +  A  phantom  incorporating  120  pm  diameter  nylon  filaments  arranged  in  rows  with  spacing 
between  fiWnts  ranging  from  0.5  mm  to  16  mm  was  constructed  and  used  to  measure  the 
point-spread-  function  (PSF)  at  various  locations  within  the  measurement  tank  as  well  as  at  a 
number  of  discrete  frequencies  within  the  envelope  of  frequencies  used  for  wide-band  data 
acquisitions.  These  measurements  were  carried  out  with  the  1  MHz  SCT  ring  in  water  and  thus 
represent  the  idealized  system  performance. 

In  general,  if  the  2-dimensional  system  point-spread  function  (the  image  of  a  point  object  or 
de  te-function  input)  is  given  by  PSF(x,y),  then  an  object  with  true  spatial  characteristics  o(x,y) 
will  be  blurred  so  that  the  image  is  represented  by  the  convolution  integral; 


f(x,y)  =  dx'  dy'o{x'ry'  )PSF{x  -x\y-y') 


which,  by  the  convolution  transform  theorem  is  equivalent  to  the  statement: 

l{u,v)  =  d{u,v)*H{u,v)  (2) 

where  I,  O  and  H  are  Fourier  transforms  of  i,  o  and  PSF  in  (1)  with  u,v  the  spatial  frequency 

variables  corresponding  to  x,y.  The  MTF  is  defined  by  the  ratio  of  the  blurred  to  the  true  object 
which  IS  given  by: 


\o{uM  _ 


(3) 


The  MTF  was  computed  from  measured  point-spread  functions  due  to  the  nylon  filaments  at 
several  locations  as  well  as  for  ten  discrete  frequencies  ranging  from  0.7  to  1.3  MHz.  A  typical 
result  of  such  measurements  is  shown  in  Figure  8a  which  suggests  that  the  cut-off  frequency  of 
the  MTF  is  somewhere  around  1  Ip/mm  which  is  roughly  on  the  order  of  the  theoretical  limit  for 
resolution  of  ^/4.  We  did  not  find  much  variation  of  the  MTF  with  frequency,  however,  and  this 
is  a  result  that  will  be  reviewed  further.  In  particular,  the  PSF  image  is  represented  by  very  few 
points  which  makes  the  MTF  calculation  (using  the  FFT)  less  accurate.  We  will  explore  the 
feasibility  of  creating  a  speed  of  sound  phantom  with  a  sharply  defined  edge  for  which  the  line- 
spread-:toction  (LSF)  can  be  measured  at  some  angle  to  provide  more  data  points  for  the  MTF 
calculation,  a  method  commonly  employed  in  x-ray  CT. 

The  effect  on  the  system  resolution  of  the  presence  of  an  aberrating  or  scattering  medium 
w^  initially  explored  by  measuring  the  image  of  a  nylon  filament  passing  through  balloons  filled 
with  solutions  of  various  speeds  of  sound.  Also  shown  in  Figure  8b  is  the  MTF  for  such  a 
measurement  and  it  can  be  seen  that  the  cut-off  frequency  has  been  noticeably  reduced  from  that 
for  the  filament  in  water.  Thus,  it  is  expected  that  the  system  MTF  in  more  severely  scattering  or 
aberrating  media  such  as  agar  and  tissue-mimicking  phantoms  will  be  degraded  even  further. 

This  topic  will  be  revisited  in  experiments  planned  for  the  0.5  MHz  SCT  ring  as  well  in  the  near 
future. 


(a) 


(b) 


0  0.5  1  15  2 

-0-  .2. 


Figure  8.  MTF  of  1  MHz  SCT  Ring.  First  Order  Bom  Magnitude  Reconstruction. 
Object  was  120  pm  dia.  Nylon  Filament  (a)  in  Water,  (b)  inside  2.75  in.  dia.  Balloon  Filled  with 
Saline  Solution  with  Speed  of  Sound  4.7%  Higher  than  Water. 
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2.4  Speed  of  Sound  Contrast  and  Linearity 

The  solutions  to  the  inverse  scattering  problem  based  on  the  Bom  approximation  result  in 
the  reconstruction  of  the  complex  propagation  constant; 

^  (^/ 3/)  ~  ^0  (^/ y)[^  +  (4) 

^^dlereko  =  (o/co;  n,  the  refractive  indexes  g™  n(x,y)=c(x,y)/co  where  c(x,y)  is  the  speed 
o  sound  of  the  object  and  cq  that  of  the  medium.  The  real  part  of  n  reconstructs  c(x,y),  the 
imaginary  part  is  an  attenuation  term  (i.e.  n(x,y)  =  a  +  ip  ). 

Measurements  of  the  linearity  of  the  reconstructed  speed  of  sound  over  a  range  in  speed 
from  - 1 1  %  to  +  9%  relative  to  that  of  water  were  carried  out  with  latex  balloon  phantoms  filled 
with  vanous  ethanol  (lower  sound  speed)  and  saline  (higher  sound  speed)  solutions.  The 
balloons  were  nominally  3/4  in.  in  diameter  thus  representing  scatterers  of  dimensions  of  several 
wavelengths  at  both  nominal  SCT  ring  frequencies  of  0.5  and  1.0  MHz.  Simulation  calculations 
by  GP  Otto  had  suggested  non-linear  behavior  of  the  reconstmcted  sound  speed  and  attenuation 
values  for  the  various  reconstmction  algorithms  currently  employed.  These  are  a  beam-forming 
algorithm  that  results  in  a  magnitude  image  (BFA)  and  two  approaches  to  corrections  to  the  first- 

order  Bom  approximation  (SM  and  STT  algorithms)  the  discussion  of  which  is  beyond  the  scope 
of  this  summary  report. 

Measurements  were  conducted  with  both  SCT  rings  by  suspending  four  solution-filled 
balloons  and  obtaining  an  image.  For  each  image,  one  of  the  balloons  was  filled  water  from  the 
scanning  tank  as  a  reference.  The  reconstmcted  images  were  analyzed  by  recording  the  contrast 
(pixel  values  relative  to  that  of  the  medium)  of  the  images  of  the  various  solutions  for  both  the 
real  and  imaginaiy  images.  An  example  of  such  an  image  pair  is  shown  in  Figure  9  which 
depicts  the  real  (9a)  and  imaginary  (9b)  parts  for  low-to-medium  contrast  solutions  and 
qualitatively  satisfies  the  notion  that  higher  sound  speed  should  be  bright  and  that  the  attenuation 
images  should  all  appear  similar  since  no  significant  differences  in  the  attenuation  values  for  the 
various  solutions  were  expected. 

Examples  of  speed  of  sound  and  attenuation  contrast  data  for  the  1  MHz  SCT  ring  are 
provided  in  Figure  10.  Similar  results  have  been  obtained  recently  with  the  0.5  MHz  ring  using 
the  same  sound  speed  solutions.  While  conclusions  are  premature  at  this  point,  it  appears  that  as 
the  reconstmction  algorithms  become  more  sophisticated,  incorporating  time-of-flight  and 
multiple  scattering  corrections,  the  linearity  range  of  both  speed  of  sound  and  attenuation  appears 
to  increase  while  contrast  or  sensitivity  decreases.  It  should  be  noted  that  “attenuation”  here  does 
not  refer  to  actual  pressure  amplitude  or  intensity  reduction  but  simply  the  image  pixel  values  for 
the  imaginary  part  of  the  image  plotted  against  phantom  sound  speed.  We  have  already  initiated 
experirnents  to  measure  similar  data  when  the  speed  of  sound  objects  are  embedded  in  strongly 
aberrating  media  such  as  graphite-loaded  agar  and  foam  or  breast-tissue  mimicking  materials  that 
we  have  been  experimenting  with.  The  intent  here,  of  course,  is  to  see  how  the  speed  of  sound 
contrast  and  linearity  are  affected  when  objects  of  dimensions  on  the  order  of  several 
wavelengths  are  embedded  in  materials  that  significantly  scatter  and  attenuate  the  incident  wave. 
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Pixel  Value  (arbitraty) 


Figure  9.  Real  and  Imaginary  Components  of  Reconstructed  Images  of  Medium  Contrast 
Balloons  using  Beam  Form  Algorithm  (TTC).  Clockwise  from  top  left;  Water,  1  73% 
Saline,  1.89%  Saline,  -1.43%  Ethanol. 
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Figure  10  a.  Speed  of  Sound  Contrast  Sensitivity  for  1  MHz  SCT  Ring  for  Three 
Reconstruction  Algorithms. 
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Figure  10  b.  Attenuation  Contrast  Sensitivity  for  1  MHz  SCT  Ring 
for  Three  Reconstruction  Algorithms. 

Examples  of  speed  of  sound  and  attenuation  contrast  data  for  the  1  MHz  SCT  ring  are 
provided  in  Figure  10.  Similar  results  have  been  obtained  recently  with  the  0.5  MHz  ring  using 
the  same  sound  speed  solutions.  While  conclusions  are  premature  at  this  point,  it  appears  that  as 
the  reconstruction  algorithms  become  more  sophisticated,  incorporating  time-of-flight  and 
multiple  scattering  corrections,  the  linearity  range  of  both  speed  of  sound  and  attenuation  appears 
to  increase  while  contrast  or  sensitivity  decreases.  It  should  be  noted  that  “attenuation”  here  does 
not  refer  to  actual  pressure  amplitude  or  intensity  reduction  but  simply  the  image  pixel  values  for 
the  imaginary  part  of  the  image  plotted  against  phantom  sound  speed.  We  have  already  initiated 
experiments  to  measure  similar  data  when  the  speed  of  sound  objects  are  embedded  in  strongly 
aberrating  media  such  as  paphite-loaded  agar  and  foam  or  breast-tissue  mimicking  materials  that 
we  have  been  experimenting  with.  The  intent  here,  of  course,  is  to  see  how  the  speed  of  sound 
contrast  and  linearity  are  affected  when  objects  of  dimensions  on  the  order  of  several 
wavelengths  are  embedded  in  materials  that  significantly  scatter  and  attenuate  the  incident  wave. 


Scatter  Amplitude  Distribution  Analysis 


Wave  propagation  in  the  breast,  according  to  a  number  of  researchers,^^'^^  is 
characterized  principally  by  scatter  and  refraction  and  these  processes  are  the  primary  causes  of 
amplitude  wavefront  distortion.  In  the  case  of  weak  scattering,  phase  distortion  can  be 
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adequately  addressed  by  various  phase-aberration  correetion  schemes.  Amplitude  distortion, 
however,  is  observed  when  scattering  becomes  strong  and  its  appearance  can  be  thought  to  define 
the  limit  of  applicability  of  weak-scattering  theories.  In  other  words,  scattering  and  refraction 
can  be  strong  and  thus  one  would  not  expect  the  first  order  Bom  approximation  to  hold  since  it 
implicitly  ^sumes  that  scattering  is  weak.  The  evidence  for  this  transition  from  weak  to  strong 
scattenng  is  the  shape  of  the  observed  probability  density  function  (pdf)  of  the  complex 
wavefront  the  modulus  of  which  is  approximated  by  a  Rayleigh  distribution.  If  scattering  is 
weak,  the  observed  distribution  should  be  Rician  and  ultimately  Gaussian. 

Given  that  the  reference  waveform  in  the  homogeneous  medium  at  the  receiver  is: 

Eoix)  =  4,(^)exp(/Vo(x))  (5) 

and  that  with  the  object  in  the  path: 

E{x)  =  A{x)q\^<P{x))  (6) 

the  wavefront  amplitude  distortion  term  is  defined  as: 

a(x)^A(x)/A^(x)  (7) 

By  forming  a  histogram  of  the  a(x)  values,  pdf  vs.  amplitude  plots  are  obtained  and 
compared  to  Rayleigh  distributions  of  identical  mean  values.  If  there  is  close  agreement  between 
the  two,  that  is  if  the  standard  deviations  of  the  computed  Rayleigh  and  the  experimentally 
determined  distributions  are  close,  then  the  scatter  field  is  considered  to  be  fully  developed. 

The  general  form  of  the  distribution  to  which  our  data  were  fitted  is  given  by  the 
expression:'^*’ 


P(A)  =  - 


r 


expi  - 


T' 


j'.(  ) 


2<t 


(3) 


where  P(A)  is  the  probability  density  for  amplitude  A,  Aq  has  integral  values  and  when  Aq  =  0, 
(8)  tends  toward  a  Rayleigh  distribution  whereas  for  large  values  of  Ao,  P(A)  becomes  Gaussim, 
lo  IS  the  modified  Bessel  function.  A  more  detailed  development  of  the  statistics  of  amplitude 
scattering  is  presented  in  Appendix  2. 

We  have  extended  this  sort  of  analysis  to  lower  frequencies  by  studying  the  scatter 
amplitude  distributions  from  several  patients  who  had  been  imaged  as  part  of  our  clinical  trials 
series  at  the  UCSD  Center  for  Women’s  Health.  Patients  had  been  imaged  initially  with  the  1 
MHz  SCT  ring  and  later  on  with  the  new  low-frequency  ring  as  well.  One  slice  from  each 
patient  was  analyzed  according  to  the  procedure  outlined  in  Appendix  2.  The  initial  analysis 
used  only  those  receivers  of  the  transducer  ring  that  were  in  the  shadow  of  the  object.  This  is 
essentially  comparable  to  the  technique  used  by  Steinberg  and  co-workers  who  used  a  single 
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transmitting  transducer  and  a  linear  receiving  array  to  record  only  pulses  transmitted  through  the 
object.  In  our  case  we  restricted  the  number  of  receiving  transducers  to  the  minimum  value 
determined  by  an  examination  of  the  sinogram  of  the  data.  We  are  also  analyzing  scatter 
amplitude  date  from  all  receivers  surrounding  the  object  and  will  compare  that  with  the  sort  of 
analysis  exemplified  in  Figure  11. 

It  is  observed  that  the  index  is  lower  for  a  breast  that  had  had  a  lumpectomy  with 
attendant  radiation  treatment  than  for  a  dense  breast  (patients  5  and  8,  respectively)  and  also  the 
index  is  lower  for  a  fatty  breast  compared  to  a  normal  breast.  Since  these  data  are  from  only  one 
coronal  slice,  it  is  not  clear  whether  this  ranking  will  hold  when  data  from  several  tomographic 
slices  of  each  breast  are  averaged.  While  the  behavior  with  frequency  for  the  low-frequency  ring 
appears  to  be  counter  to  expectations,  the  classifying  index  based  on  the  mean  value  and  standard 
deviation  of  the  recorded  amplitude  distribution  (i.e.  p  ==  <x>/a )  appears  to  differentiate 
between  breast  tissue  types. 
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Figure  1 1.  Scatter  Distribution  Index  |j,  —  <x>/or  as  a  Function  of  Frequency  for  Four 
Patients  with  Breast  Types:  Pat.  5:  Dense;  Pat.  8:  Fatty/Rad.Tx;  Pat.  14:Normal;  Pat.  15:Fatty 

These  analyses  and  further  results  are  described  in  more  detail  in  the  paper  of  Appendix  3. 
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2.6  Tissue  Contrast  Determination 


Ultrasound  diffraction  tomography  is  considered  a  quantitative  imaging  technique  since 
there  exists  a  direct,  albeit  complex,  relationship  between  the  nature  of  the  scatter  field 
(amplitude  and  phase)  and  the  magnitude  of  the  fluctuations  in  tissue  properties  such  as 
compressibility,  density  and  attenuation.  In  principle,  such  direct  relationships  hold  for  the  case 
o  exact  solutions  to  the  wave  equation.  It  is  not  clear  over  what  magnitude  ranges  in  parameter 
actuations  the  images  are  quantitative  when  linearizing  approximations  are  used  to  solve  the 
wave  equafion  and  when  noise,  artifacts  and  aberrations  are  present.  To  make  an  initial  empirical 
etermmation  of  the  nature  and  magnitude  range  of  image  parameters,  a  number  of  patient 
images  from  our  clinical  trials  were  examined  in  detail  in  order  to  ascertain  the  following: 


•  Determination  of  tissue  discrimination  ability  by  checking  values  and  consistency  of  image 
contrast  of  different  tissue  types  within  the  breast  and  between  patients. 

Determination  of  the  uniformity  of  the  pixel  or  gray  scale  value  for  key  image  features  (e  g 

water )  within  the  breast,  for  different  patients,  the  two  transducer  rings  and  various 
reconstruction  algorithms. 

presentations  (real,  imaginary,  magnitude)  and  reconstructions  (BFA, 
oMA,  b  1 T)  which  yield  optimal  tissue  discrimination  and  consistency. 


2.6.1,  Approach 


This  is  an  initial  attempt  at  establishing  a  procedure  for  analyzing  images  and  as  such  it  is 
subject  to  change  and,  hopefully,  improvements.  The  evaluation  procedures  discussed  here  deal 
not  with  the  raw,  time-dependent  data  but  with  the  reconstructed  images  and  the  pixel  values  of 
image  elements  of  clinical  significance.  A  radiologist  (Dr.  Ysrael  in  most  instances)  viewed 
images  of  coronal  breast  slices  on  the  computer  monitor  and,  with  the  aid  of  corresponding 
mammograms,  identified  the  locations  (x/y  coordinates)  and  extents  of  significant  tissue 
structures  such  as  :  Skin,  Subcutaneous  Fat,  Parenchyma,  Ligaments,  Cysts,  Lesions  and 
Unidentified  Features.  Examples  of  the  last  item  would  be  bright  dots  or  specks  that  might 
represent  calcifications  or  some  unidentified  strongly  reflecting  (or  scattering)  small  object. 

0  owing  this  identification  of  key  image  features  each  was  evaluated  by  recording  mean  pixel 
values  for  regions  of  interest  (ROTs)  that  encompassed  the  feature.  For  each  image  slice,  the 
mean  pixel  value  for  water  was  obtained  by  defining  an  annular  region  that  excluded  the  breast 
object  as  well  as  artifact-filled  regions  adjacent  to  the  transducer  ring.  The  gray-scale  value  for 
w^er  was  used  to  compute  the  contrast  for  each  image  element  from  the  relationship: 
c  ^object  -  Pwater )/  Ipwaterl-  The  use  of  the  absolute  value  for  the  water  pixel  value  assures  that 
t  e  sign  of  the  contrast  value  unambiguously  indicates  the  relationship  to  the  pixel  value  of  water 
which,  very  often,  happens  to  be  less  than  zero. 

2.6.2.  Preliminary  Results 

Image  slices  from  eight  patients  (5,7,9,12,13,14,19  and  21)  obtained  with  the  two 
transducer  nngs  were  analyzed  in  this  fashion.  Since  generally  10  or  so  tomographic  slices  are 
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available  for  each  patient  and  since  some  were  recorded  with  both  SCT  rings,  even  this 
preliminary  data  set  is  fairly  large  and  time-consuming  to  produce.  From  a  recording  of  the 
dimensions  of  the  transducer  ring,  the  pixel  dimension  is  calculated  and  used  to  compute  the 
sizes  of  objects  or  features  identified  by  the  radiologist.  Thus,  for  the  1  MHz  ring  (Patients  5  7  9 
and  12),  pixel  dimension  are  0.21  x  0.21  mm  and  for  the  fifth  slice  of  Patient  5,  for  example’  the 
overall  breast  dimensions  are  98.25  x  99.30  mm.  Although  the  pixel  value  of  the  skin  layer  was 
recorded  for  many  cases,  it  was  found  to  range  quite  a  bit  and  is  not  believed  to  be  of  diagnostic 
value.  It  is  principally  a  function  of  the  measurement  geometry  and  may  be  an  as  yet  poorly 
understood  edge  artifact. 

a.  Infra  Breast  Data 

One  of  the  important  questions  that  was  to  be  answered  by  this  sort  of  review  of  the  data 
is  the  consistency  of  imap  elements  from  slice  to  slice.  Clearly,  if  the  same  feature  or  tissue 
type  vanes  significantly  in  contrast  from  slice  to  slice,  there  is  little  hope  of  ever  achieving 
general  quantitative  guidelines  that  will  permit  tissue  identification  on  the  basis  of  relative 
contrast  level  (i.e.  scattering  strength).  In  order  to  asses  this,  all  slices  for  patients  5 , 7,  13  and 
14  were  analyzed  and  evaluated.  The  results  are  presented  in  the  attached  charts  and  some  of  the 
data  are  discussed  briefly  to  provide  insight  into  the  nature  of  the  contrast  data. 

Patient  5  is  a  dense,  normal,  right  breast  with  diffuse  calcifications.  The  plot  of  the  raw 
pixel  value  data  points  reveals  values  much  higher  than  for  the  other  patients  as  well  as  some 
variation  in  the  water  value  from  slice  to  slice.  The  general  trend  is  for  the  water  values  to 
decrease  toward  the  chest  wall  and  that  may  in  part  be  due  to  the  change  in  the  area  of  the 
annular  ring  over  which  water  values  are  averaged.  The  contrast  values  for  parenchyma  tissue  are 
around  40%  higher  than  water  behind  the  nipple  and  then  jump  to  around  100%  above  water  in 
the  middle  of  the  breast.  Subcutaneous  fat  generally  is  lower  than  water  by  as  much  as  50  %. 
These  trends  are  reflected  in  the  mean  object  contrast  which  is  reasonably  low  for  the  first  three 
slices  but  then  jumps  to  values  of  40  %  or  greater  in  the  center  of  the  breast.  Figure  12  is  a  plot 
of  percent  tissue  contrast  for  Patient  5  over  the  whole  volume  of  the  breast. 

We  observe  in  general  that  while  differences  between  contrast  values  of  different  tissue 
types  are  fairly  consistent  throughout  the  breast  volume,  there  appears  to  be  a  trend  for  all 
contrast  values  to  increase  toward  the  chest  wall.  It  may  be  that  this  is  somehow  a  result  of  the 
fact  that  the  object  dimensions  increase  toward  the  chest  wall  and  thus  the  total  attenuation 
through  the  object  increases  as  well.  Phantom  studies  will  have  to  be  carried  out  to  determine  the 
exact  nature  of  this  behavior. 

The  contrast  values  for  subcutaneous  fat  are  generally  at  or  below  the  level  of  water 
suggesting  that  for  medium-matching  experiments  one  may  wish  to  lower  the  water  speed  of 
sound  with  alcohol.  However,  how  one  changes  the  medium  sound  speed  may  well  depend  on 
what  one  wants  to  emphasize.  Imaging  the  subcutaneous  fat  may  be  not  nearly  as  important  as 
imaging  the  glandular  tissue  and  structures  contained  within  it.  Thus,  one  may  want  to  use  saline 
solutions  to  increase  the  medium  speed  of  sound  to  more  closely  match  the  overall  object  value. 
Indeed,  this  is  an  interesting  point  for  further  investigation,  for  if  one  were  to  match  the  fat, 
penetration  through  that  layer  would  be  enhanced  but  then  the  overall  contrast  of  the  glandular 
tissue  would  even  be  higher  and  thus  any  advantage  gained  may  be  lost  in  tbe  reconstruction. 
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PAT  5  CONTRAST  DATA 


Figure  12.  Percent  Tissue  Contrast  for  Normal  Breast  Tissue  Elements  for  Patient  5. 


All  0.5  MHz  SCT  ring  images  analyzed  were  first-order  time-of-flight  corrected  images 
(BFG)  and  have  real  and  imaginaiy  components.  The  overall  object  contrast  values  for  some  of 
the  patients  are  quite  high,  exceeding  100%.  This  in  itself  may  be  an  indicator  that  the  breast  is 
either  generally  quite  dense  or  contains  significant  elements  of  high  density. 

b.  Inter-Breast  Comparison 

The  relative  consistency  of  breast  tissue  contrast  within  the  breast  is  a  critical  requirement 
but  also  clearly  desirable  would  be  some  rough  correspondence  between  the  same  tissue  type  for 
different  patients.  To  the  extent  that  this  was  found  to  be  the  case  for  the  patients  and  particular 
reconstructions  evaluated  is  shown  in  Figure  13.  Here  are  presented  the  magnitude  tissue 
contrast  values  for  five  patients  (5,7,9,12,13)  and  this  comparison  is  of  considerable  interest. 

The  data  in  Figure  13  are  for  normal  tissue  components  such  as  fat  and  parenchyma  and 
two  general  observations  can  be  made:  The  trend  for  contrast  values  to  increase  toward  the  chest 
wall  is  reflected  in  all  cases  and  although  the  data  exhibit  considerable  scatter,  groupings  by 
tissue  type  for  all  patients  are  quite  clearly  discernible.  This  suggests  then,  on  the  basis  of  veiy 
limited  data,  that  the  ultrasound  diffraction  images  are,  indeed,  quantitative  and  able  to 
distinguish  between  normal  breast  tissue  types. 
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1  MHZ  DIFF.  CONTRAST  SUMMARY 
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Figure  13.  Summary  of  1  Mhz  SCT  Data  of  Normal  Breast  Tissue  Contrast  for  5 
Patients.  Note  that  Parenchyma  is  open  symbols,  Fat  shaded  symbols  and  the 
mean  Object  contrast  is  represented  by  solid  symbols. 


A  summary  comparison  of  data  from  eight  patients  is  afforded  by  Figure  14  which 
groups  both  normal  and  pathological  tissue  types  whose  values  have  been  averaged  over  the 
whole  breast  volume.  This  data  were  obtained  with  both  the  0.5  MHz  and  the  1  MHz  SCT  rings. 
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8  PATIENT  DATA  SUMMARY 
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Figure  14.  Summary  of  Mean  Breast  Tissue  Contrast  Values  for  Eight  Patients  Obtained 
with  Two  Transducer  Rings  (0.5  MHz  and  1  MHz). 


bl.  Object  Contrast 

The  overall  object  contrast  is  seen  to  vary  quite  a  bit  and  is  surprisingly  low  for  Patient  7, 
something  that  a  cursory  evaluation  of  the  images  would  not  suggest.  It  is  to  be  noted  that  the 
values  presented  here  are  means  for  all  slices  and  these  can  fluctuate  quite  a  bit  depending  on  the 
distribution  of  fat  vs.  parenchyma,  for  example.  As  we  saw  in  the  previous  section,  there  appears 
to  be  a  tendency  for  the  contrast  values  of  most  breast  tissue  components  to  be  highest  toward  the 
center  of  the  breast. 

With  the  possible  exception  of  Patient  7,  the  overall  object  contrast  values  are  so  high  that 
serious  concerns  about  the  applicability  of  the  Born  Approximation  have  to  be  raised.  With  such 
high  speed  of  sound  variations  and  large  object  sizes  total  phase  shifts  with  the  1  MHz  ring  can 
readily  exceed  15  wavelengths  and  would  be  half  that  for  the  lower  frequency  ring.  Since  these 
are  magnitude  images  it  is  difficult  to  determine  what  component  of  contrast  is  due  to  speed  of 
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speed  of  sound  and  what  is  due  to  attenuation  (we  assume  constant  density  since  it  is  not 
explicitly  contained  in  the  wave  equation).  This  suggests  that  the  diffraction  images  are  not 
explicitly  speed  of  sound  nor  attenuation  images  but  rather,  images  of  fluctuations  of  the 
propagation  constant,  k.  That  is,  k  is  a  complex  number  whose  real  component  is  the  speed  of 
sound  and  whose  imaginary  component  is  the  attenuation.  With  large  objects  or ,  more 
specifically,  with  large  phase  shifts,  these  components  become  mixed  as  a  result  of  phase 
wrapping. 

hi.  Parenchyma 

The  contrast  value  for  parenchyma  is  consistently  high  and  quite  consistent  for  the  two 
rings.  That  it  seems  about  twice  as  high  with  the  .5  MHz  ring  as  with  the  1  MHz  ring  may  well 
be  a  reflection  of  frequency  and  not  due  to  significant  tissue  differences  between  patients.  Such 
questions,  of  course,  are  most  easily  addressed  by  scanning  the  same  breast  with  both  rings, 
everything  else  being  as  equal  as  possible.  We  have  attempted  this  with  Patient  22/23  but  system 
problems  have  prevented  us  from  reaching  any  conclusions. 

b3.  Subcutaneous  Fat 

The  tendency  for  subcutaneous  fat  to  have  lower  contrast  than  water  is  consistent  for  all 
patients  and  again,  we  observe  that  the  lower  frequency  transducer  ring  produces  more  than  twice 
the  contrast  value  observed  at  higher  frequency. 


b4.  Cyst 

Although  only  one  data  point  exists  so  far  (Patient  7),  it  is  encouraging  to  note  that  the 
contrast  value  for  a  large  fluid-filled  cyst  is  close  to  that  of  water,  something  that  was  not 
apparent  from  a  visual  examination  of  the  image.  Additionally,  septations  suggested  by 
ultrasound  B-scan  images  were  clearly  identifiable  in  the  diffraction  image. 

b5.  Cancer 

Contrast  values  for  cancer  appear  quite  high,  somewhere  between  30  and  60  percent 
above  that  of  parenchyma.  In  addition  to  the  higher  numerical  value,  the  appearance  of 

p^enchyma,  when  a  large  enough  region  of  interest  is  selected,  has  distinct  attributes  that  may  be 
diagnostic. 
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2.7  Publications  and  Presentations 


Some  of  the  work  of  this  fellowship  was  aecepted  for  publieation  or  presentation  by  the 
journals  and  conferences  listed  below.^^"*’  This  effort  was  in  all  cases  the  collaboration  of 
several  individuals  and  institutions  and  the  opportunity  to  participate  in  these  efforts  is  gratefully 
acknowledged. 

•  HS  Janee,  MP  Andre,  PJ  Martin,  GP  Otto,  BA  Spivey,  JP  Jones,  "Performance 
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2 1(6):  1002,  1994. 
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CONCLUSIONS 


Although  this  is  the  final  report  of  the  pre-doctoral  fellowship,  the  research  initiated 
within  this  program  is  ongoing  and  many  conclusions  are,  perforce,  preliminary.  This  is  not  a 
surprising  situation  for  a  new  research  area  that  explores  a  unique  approach  to  breast  imaging  and 
requires  not  only  new  interpretations  of  images  but  an  adaptation  to  a  new  presentation  format. 
Some  aspects  of  the  fellowship  grant,  in  particular  those  dealing  with  training  and  certification, 
require  no  conclusions  other  than  the  statement  that  certain  objectives  were  realized  and  others 
are  left  to  be  accomplished. 


Certification  in  Diagnostic  Radiological  Physics/Mammography 


A  considerable  body  of  experience  in  diagnostic  imaging  equipment  quality  assurance 
and  control  has  been  obtained  through  course  work  at  the  University  of  California,  Irvine  and 
through  surveys  of  many  systems  at  a  number  of  locations  under  the  supervision  of  a  board- 
certified  medical  physicist.  Special  efforts  were  made  to  prepare  for  certification  in 
mammography  with  the  State  of  California  and  an  application  for  certification  has  been 
submitted.  With  the  granting  of  the  M.S.  degree  in  Radiological  Sciences,  experience  gained 
under  the  supervision  of  a  certified  physicist  can  be  counted  toward  the  work-experience 
requirements  for  eligibility  to  take  the  qualifying  board  examinations  offered  by  the  AAPM  or 


II.  Dissertation  Research 


With  the  addition  of  the  low-frequency  transducer  ring,  many  of  the  performance 
charactenzation  tasks  had  to  be  carried  out  for  both  units.  Thus,  the  Statement  of  Work  has  been 
expanded  in  scope  if  not  in  the  nature  of  the  tasks  proposed.  Ultimately  the  decision  may  be 
made  to  rely  on  only  one  of  the  transducer  rings  for  general  breast  imaging.  At  this  point,  not 
enough  information  about  the  relative  advantages  and  drawbacks  of  the  two  rings  is  available  yet 
to  formulate  such  a  decision. 

The  measurements  of  transducer  radiation  patterns  and  tomographic  slice  profiles  carried 
out  so  far  suggest  that  the  transducers  of  the  1  MHz  SCT  ring  produce  a  radiation  pattern  that  is 
a  fair  approximation  to  a  dipole  pattern.  This  cannot  be  said  for  the  lower  frequency  ring  the 
transducers  of  which  have  a  lower  L/D  ratio  and  thus  are  characterized  by  a  more  complex 
radiation  pattern.  It  will  be  one  of  the  tasks  for  next  year  to  investigate  the  effect  of  the  deviation 
from  a  dipole  pattern  on  the  reconstructed  images.  It  is  to  be  noted  here  that  the  measurements 
so  far  have  looked  only  at  a  few  individual  transducer  elements  and  that  it  is  assumed  that  the 
receive  pattern  of  the  transducers  is  similar  to  their  transmit  pattern.  There  are  indications  from 
the  examination  of  the  "noise  image",  that  is  the  image  of  the  water-filled  tank,  that  these 
assumptions  are  largely  valid. 
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The  effective  slice  thickness  has  been  found  to  be  close  to  the  long  dimension  of  the 
transducers,  an  observation  confirmed  by  examining  the  image  of  a  nylon  filament  ramp  obtained 
with  the  1  MHz  SCT  ring.  For  clinical  applications  to  breast  imaging  that  slice  width  is 
undesirably  wide  and  considerations  are  being  given  to  perhaps  a  combination  of  hardware  as 

well  as  data  acquisition  and  processing  procedural  changes  that  might  reduce  the  effective  slice 
width  of  the  image. 

System  MTF  measurements  to  date  have  been  only  carried  out  for  the  higher  frequency 
transducer  ring  and  suggest  that  the  cut-off  frequency  is  in  reasonable  agreement  with  the 
theoretically  expected  resolution-limit.  This,  however,  holds  only  for  the  ideal  case  for  which  a 
point-scatterer  is  suspended  in  the  reference  medium(water).  The  MTF  becomes  degraded  when 
an  object  such  as  a  cylinder  filled  with  a  somewhat  different  sound-speed  solution  is  introduced 
and  thus  more  serious  deterioration  is  expected  when  a  strongly  scattering,  large  object  like  a 
breast  is  imaged.  This  effect  will  be  explored  in  future  experiments  as  will  be  the  performance  of 
the  0.5  MHz  transducer  ring.  The  investigation  of  the  changes  in  MTF  with  transducer  frequency 
so  far  has  proven  inconclusive.  Some  cyclical  variations  with  frequency  are  observed  that  are  not 
yet  understood  but  may  reflect  a  basic  system  operation/image  reconstruction  artifact. 

Speed  of  sound  contrast  and  linearity  measurements  have  been  carried  out  for  both 
transducer  rings  but  only  the  data  for  the  1  MHz  ring  have  been  analyzed  so  far.  These 
measurements  were  for  homogeneous,  cylindrical  objects  that  were  several  wavelengths  in 
diameter  and  ranged  in  speed  of  sound  relative  to  that  of  water  from  around  -11%  to  +8%.  The 
contrast  of  the  reconstructed  images  when  plotted  against  the  refractive  index  or  change  in  sound 
speed  has  linear  ranges  the  widths  of  which  (in  terns  of  sound  speed)  depend  on  the 
reconstruction  algorithm  employed.  The  simplest,  first-order  Bom  approximation  reconstmction 
yields  the  lowest  linear  range  with,  however,  the  highest  sensitivity  (defined  as  the  change  in 
image  contrast  with  sound  speed).  More  sophisticated  reconstructions  that  employ  time-of-flight 
corrections  and  attempts  to  treat  multiple  scattering  appear  to  have  wider  linear  ranges  but 
markedly  reduced  sensitivity.  It  is  premature  to  suggest  that  these  are  general  conclusions  and 
for  the  moment  they  only  apply  to  a  very  narrow  set  of  test  objects  and  experimental  conditions. 

In  a  departure  from  the  Statement  of  Work  a  task  was  added  that  examined  the  scattering 
process  in  breast  tissue  in  some  detail.  A  statistical  analysis  of  the  raw  scattering  data  revealed 
that  the  amplitude  distributions  of  these  data  appear  to  correlate  with  breast  tissue  type,  though 
the  exact  nature  of  that  correlation  has  not  yet  been  established.  If  this  trend  holds  reliably  for 
many  patients  and  breast  types  then  it  may  be  possible  to  design  data  acquisition  and  image 
reconstmction  procedures  that  are  optimized  for  each  patient  based  on  a  quick  statistical  analysis 
of  the  scattered  data  from  a  precursor"  insonification  of  the  tissue.  In  a  preliminary  analytical 
effort,  an  operator-based  formalism  has  been  developed  that  suggests  that  a  range  of  linearizing 
approximations  to  the  wave  equation  are  possible  and  that  the  Bom  and  Rytov  approximations 
represent  opposite  endpoints  to  this  range.  Although  the  implementation  of  reconstmctions 
based  on  such  a  general  formalism  are  not  obvious  and  beyond  the  scope  of  this  work,  the 
preliminary  results  of  the  statistical  scatter  analysis  suggest  that  different  breast  images  may  be 
optimally  reconstructed  by  a  more  Bom-like  or  more  Rytov-like  approximation.  This  is  an 
exciting  development  the  difficulty  of  implementation  of  which  cannot,  however,  be 
underestimated. 

A  second  modification  of  the  Statement  of  Work  addresses  the  quantitative  assessment  of 
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tissue  contrast  for  both  normal  breast  tissue  and  pathologic  tissue.  Analogous  to  the  Hounsfield 
scale  used  in  x-ray  computed  tomography,  tissue  contrast  is  referenced  to  water  and  appears  to  be 
reliably  differentiated  by  tissue  type  such  as  fat,  parenchyma  ligaments  and  vessels.  A  very 
encouraging  preliminary  finding  is  that  benign  structures  such  as  fibrocysts  and  fibroadenomas 
and  malignant  lesions  can  be  quantitatively  differentiated  from  normal  breast  tissue.  It  is, 
however,  critically  important  to  stress  that  these  findings  are  based  on  a  very  small  sampling  of 
the  population  and  general  conclusions  must  be  deferred  until  a  mueh  larger  and  comprehensive 
data  base  is  available. 
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APPENDIX  1 


DIAGNOSTIC  IMAGING  QUALITY  CONTROL  EXPERIENCE  SUMMARY 


DATE 

— 

■QC  ACTIVITY 

8-7-92 

VAMC* 

GE  X-ray  CT  Unit  Annual  Survey 

9-8-93 

ft 

GE  Senographe  DMR  Mammography  Unit 
Acceptance  Testing 

6-23-94 

Children's  Hospital** 

Mobile  X-Ray  Unit  Annual  Survey 

7-23-94 

ff 

Bi-Lateral  Cardiology  Unit  (Toshiba) 
Acceptance  Testing 

10-8-94 

ff 

Mobile  X-Ray  Units(5)  Annual  Survey 

12-4-94 

VAMC 

Accreditation  Testing,  GE  DMR 

Mammography  Unit 

12-15-94 

ff 

Cardiac  Cath  Lab  Radiation  Shielding 
Evaluation 

12-16-94 

Children's  Hospital 

Radiography/Fluoroscopy  Annual  Surveys, 
Rooms  23  and  5 

12-22-94 

ff 

X-Ray  CT  Scanner,  Room  4,5  Fluoroscopy 

Units  Annual  Survey 

12-27-94 

VAMC 

GE  DMR  Mammography  unit  QC 

1-2-95 

Children's  Hospital 

OR  C-Arm  Unit,  2  Mobile  C-Arm  units. 

Mobile  Dental  Unit  Annual  QC 

3-10-95 

VAMC  Outpatient 
Center*** 

Radiography/Fluoroscopy  Suite  Shielding 
Evaluation/Survey 

5-10-95 

VAMC 

Room  3  Fluoroscopy  Unit 

Annual  QC 

5-17-95 

ff 

Room  1  Radiography  Unit 

Annual  QC 

5-25-95 

ff 

Room  1  Fluoroscopy  Unit 

Annual  QC 

6-15-95 

ff 

Room  7,  Chest  Unit 

Annual  QC 
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6-26-95 

ft 

X-Ray  CT  Scanner  Annual  QC 

6-27-95 

ff 

Room  4,  Cath.  Lab  Fluoroscopy  Unit 

Annual  CT 

1-7-96 

GE  Senographe  DMR  Mammography  Unit 

Annual  Survey 

*  Veterans  Affairs  Medical  Center,  La  Jolla,  CA. 

**  Sharpe/Children's  Hospital,  San  Diego,  CA 
***  Veterans  Affairs  Medical  Center  Outpatient  Clinic,  San  Diego,  CA. 
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APPENDIX  2 


SCATTER  AMPLITUDE  DISTRIBUTION  ANALYSIS 


A2.1  Statistics  of  the  Complex  Amplitude 


When  a  wave  sueh  as  an  ultrasonic  plane  wave  interacts  with  an  object  that  is 
inhomogeneous  and  composed  of  an  assembly  of  random  scatterers,  the  resultant  field  is 
described  by  a  complex  scattering  amplitude^^’^* ; 


=  Re{^}  =  ^  ^lajcos^^t 
=  Im{.4}  =  ^  J]|ajsin<;>i 


(A2-1) 


It  is  to  be  noted  that  the  average  values  of  these  amplitudes  over  a  macroscopic  ensemble  is  0. 
The  real  and  imaginary  parts  of  the  complex  field  thus  have  zero  means,  have  identical  variances 
and  are  uncorrelated. 

If  the  number  N  is  very  large,  i.e.if  there  are,  as  in  our  application,  many  scatterers,  and 
the  field  thus  is  the  sum  of  many  independent  random  variables,  then  by  the  Central  Limit 
Theorem  ’  ,  the  amplitude  distribution  is  Gaussian: 

The  pressure  scatter  field  p  is  given  by: 


p  =  Aexp(i^)  =  X  +  zY 


(A2-2) 


where  A  is  the  amplitude,  <()  the  phase  and  X  and  Y  are  the  sums  of  many  independent 
contributions: 


X  = 


and 


so  that  the  field  is  the  sum  of: 


N  N 

p  =  ^  A.  exp(z>,)  =  ^(X,  +  iX) 


(A2-3) 


(A2-4) 
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The  Central  Limit  Theorem  says  that  the  distributon  of  a  random  variable  which  is  the  sum  of  N 
independent  random  variables  approaches  normal  as  N  becomes  veiy  large  regardless  of  the 
distribution  of  each  random  variable.  Thus  X  and  Y  are  normally  distributed  if  N  is  very  large. 
Next,  it  is  assumed  that  the  phase  <]);  is  uniformly  distributed  over  2%  since  each  component 
A,  exp(z^;)  is  random  and  may  be  assumed  to  be  uniformly  distributed  over  2n. 


The  probability  density  function  of  the  amplitude  A  and  the  phase  ^  is  then  as  follows: 

P{A,<p)=  P{A)P{^),  P{(^=\llii,  Q<(f><2n  (A2-5) 

It  can  be  easily  shown  that: 

00  2n 

(X)  =  (A  cos  =  JdA  J*d^Acos^P(A,^)  =  0 

0  0 

(y)  =(Asin^zi)  =0 

00 

(XY)  =  (a"  sin  ^cos  <()  =  0,  (x')  =  (  Y')  =  <j  =  J;rA'P(A)dA  (A2-6) 

0 

It  can  be  shown  that,  since  X  and  Y  are  normally  distributed: 

P(X,  Y)  =  P(X)P(Y)  =  ^  exp[-{X’  +  y“)/2cr"]  (A2-7) 

P{X,Y)dXdY  =  P{A,^)dAd^ 

2k 

Now  dXdY  =  AdAd<t),  P(A)  =  (A2-8) 

0 

Then:  P(A)  =  -4-exp(-A^  /  2o-")  (A2-9) 

G 

This  is  the  Rayleigh  distribution  which  is  characterized  by  a  single  parameter,  the  standard 
deviation.  Thus  the  amplitude  and  phase  of  the  scattered  waves  are  independent,  with  the 
amplitude  Rayleigh  distributed  and  the  phase  uniformly  distributed. 

The  development  of  the  Rician  or  Rice-Nakagami^^  distribution  follows  along  similar  lines. 
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Let  <p>  =  Aoexp(i(|)o)  be  the  coherent  field.  Letting  the  reference  phase  be  0,  the  incoherent 
field  is  given  by  Ps  =  p  -  <p>  and  is,  as  above,  written  as: 

N  N 

P.  =  (^-^o)  +  *X  X-Ao  =  ^X,,  ^  =  (A2-10) 

i  i 

The  probablity  density  function  W(A),  based  on  the  Central  limit  theorem,  is  now  given  by: 


PV(X,Y)  = 


1  .  r 


(x-A,)Hy^1 

J 


(A2-11) 


Then,  by  a  similar  argument  to  (A2-7  through  A2-9),  the  probability  density  for  the  amplitude  A 
is  given  by: 


W(A)=  IW(X,Yyirf«i  =  4exp(-^^^)7„(^)  (A2-I2) 

J  (7  1<J  <T 

0 

Where  Iq  is  the  modified  Bessel  function.  When  the  constant  phasor  Ao  — >0,  this  reduces  to  the 
Rayleigh  distribution  in  ( A2-9).  For  large  values  of  Ao ,  W(A)  becomes  Gaussian.  The  ratio  of 
the  mean  value  of  the  distribution  W(A)  to  its  standard  deviation  is  found  to  increase  with  Ao  and 
thus  provides  a  useful  index  to  indicate  whether  a  scatterer  is  mostly  Rayleighian  or  Rician.  The 
Tissue  Index  is  defined  as  :  TI  =  a  /cr^  and  can  be  shown  to  be  given"*^  by: 


(1  + 


—)4  (—)+—/,(—) 
2  2  2  2 


(A2-13) 


The  ratio  is  unity  for  an  exponential  distribution,  1 .24  for  a  Rayleigh  distribution  and  becomes 
arbitrarily  large  with  Ao  for  a  Rician  distribution. 


A2.2  Breast  Tissue  Scatter  Amplitude 

The  wavefront  amplitude  distortion  term  defined  in  the  Body  section  (equation  7)  for  our 
particular  experimental  geometry  is  given  by: 


1 


N-\ 


k=\ 


j^k 


a(x)  =  A(x)/AM) 

Uj  =  A^j  /  AJ  where  <  Aj>- 
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(A2-14) 


» 


where  the  index  j  refers  to  the  receivers  in  the  transducer  ring  and  k  to  the  transmitters.  As 
defined,  the  amplitude  ratio  aj  for  a  particular  receiver  location  is  thus  the  mean  of  N-1  views  of 
the  object  where  N  is  either  1024  or  512. 

The  histogram  data  of  the  ratio  of  scattered  pressure  amplitude  to  that  observed  with  the 
reference  medium  (water)  are  analyzed  to  see  whether  they  are  more  nearly  fit  with  a  Rayleigh  or 
Rician  distribution.  The  Rayleigh  distribution  in  scatter  amplitude  results  when  scattering  and 
refraction  are  strong,  suggesting  that  the  object  is  composed  of  a  large  number  of  independent, 
random  scatterers  which  produce  uniform  phase  distributions  and  normal  amplitude  distributions. 
The  modulus  of  such  a  distribution  is  Rayleighian.  As  scattering  and  refraction  become  weaker, 
the  distribution  tends  toward  a  Rician  distribution  and  ultimately  toward  a  Gaussian. 

A  Tissue  index  was  defined  as  :  m  =  <x>  /  <7  where  <x>  is  the  mean  of  the  distribution  and 
CT  its  standard  deviation.  This  ratio  is  relatively  low  for  a  Rayleigh  distribution  and  increases  as 
the  distribution  tends  toward  Rician. 

An  example  of  a  simple  analysis  using  the  Mathsoft  Mathcad  application  is  enclosed.  The 
example  is  for  a  tomographic  slice  of  the  breast  of  Patient  15  which  was  classified  by  the 
mammography  radiologist  as  dense.  It  is  seen  that  the  scatter  amplitude  distribution  is  best 
approximated  by  a  Rician  distribution  with  k  =  1 . 


40 


STATISTICAL  EVALUATION  OF  SCATTER  DATA  FROM  BREAST 
TISSUE 

PATIENT  15  -  1.0  MHZ  SCT  RING 

FILE:  Patl5d2 

a.  Rayleigh  Distribution  Fit  to  Data 
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y 


m 


m 
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Normalized  Histogram 


Determination  of  Std.  Deviation 
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Frequency 


p 


Patient  15,  f=  j  Band.  Rayleigh  Fit 


m 

Rel.  Amplitude 


Figure  1.  Scatter  Amplitude  Distribution  for  Patient  15,  Slice  6  at  nominally  0.70 
MHz.  Rayleigh  Fit  to  Distribution.  Std.  Deviation  calculated  for  normalized 
measured  distribution.  DISTRIB2  with  banding  and  with  (8:1)  smoothing. 


b.  Rician  Distribution 
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Frequency 


Patient  15,  f=  j  Band:  Rician  Distrib. 


Figure  2.  Rician  Distributions  Fit  to  Patient  15,  Slice  6  Data  from  Figure  1. 
For  Several  values  of  k.  k=0  tends  toward  Rayleigh  Distribution,  DISTRIB2 


